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l.O  INTRODUCTION 
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During  the  past  four  years  an  ONR-sponsored  research  program  has 
been  underway  at  North  Carolina  State  University  which  basically 
involves  the  use  of  the  Monte  Carlo  method  to  study  high-field 
transport  properties  of  semiconductor  materials.  The  emphasis  has 
been  directed  toward  the  study  of  ternary  and  quaternary  Ill-V  compounds 
in  an  attempt  to  identify  specific  materials  with  desirable  electronic 
properties  for  microwave  device  applications. 

A substantial  effort  has  been  devoted  to  developing  a complete 
Monte  Carlo  computer  program  which  can  be  used  as  a predictive  tool 
in  the  study  of  hot  electron  transport,  as  well  as  low  field  ohmic 
transport,  in  a general  class  of  semiconductor  materials.  We  presently 
believe  that  the  capabilities  which  have  been  developed  so  far  under 
ONR  sponsorship  are  second  to  none  within  the  scientific  community, 
and  that  the  bulk  static  transport  properties  of  two  general  classes 
of  quaternary  III-V  semiconductors  can  be  predicted  using  our  computational 

methods.  These  two  classes  of  quaternary  materials  consist  of  compounds 

IT  u r . .III„III  II1,^V  , ,III_V„V_V, 

of  the  form  A,  B C,  D and  A B C D (or  conver'-  y A BCD), 
1-x  X 1-y  y xyz  xyz 

which  can  be  specified  to  include  any  binary  or  ternary  compound  as  a 
subset  of  these  general  material  forms. 

The  accuracy  of  the  computational  program  depends  on  several  factors. 
One  major  factor  is  the  way  in  which  the  quaternary  material  parameters 
are  calculated.  The  technique  depends  on  a knowledge  of  the  binary 
material  parameters  and  either  theoretical  models  or  experimental  data 


1? 


for  the  ternary  material  parameters.  This  Information  Is  coupled  Into 
a general  Interpolation  procedure  for  calculating  the  quaternary 
material  parameters,  from  which  the  transport  properties  are  obtained 
In  the  Monte  Carlo  method.  This  work  has  culminated  In  the  prediction 
that  specific  compositions  of  the  quaternary  compound  Ga^^  ^In^Pj^_^As^ 
lattice-matched  to  InP  substrates  possess  transport  properties  having 
distinct  advantages  over  CaAs  [1]  for  device  applications.  In  the  past 
eighteen  months  considerable  experimental  work  on  this  compound  has 
shown  the  material  to  indeed  be  promising  for  many  device  applications. 
However,  the  complete  advantages  due  to  the  predicted  superior  trans- 
port properties  have  not  been  realized  at  the  present  time,  and  efforts 
to  make  a microwave  MESFET  have  only  been  partially  successful  [2). 
However,  It  appears  that  many  of  the  device  problems  can  be  overcome 
through  the  refinement  of  the  materials  technology,  and  there  are 
reasons  to  remain  optimistic  that  the  predicted  material  advantages  will 
be  borne  out.  In  the  authors'  opinion  the  influence  of  alloy  scattering 
still  remains  an  unanswered  question.  However,  very  recent  results  on 
AI^  ^Ga^As  strongly  suggest  that  our  approach  which  uses  the  electron 
affinity  difference  for  the  alloy  scattering  potential  is  correct  [3). 

It  Is  to  be  noted  that  these  recent  results  have  come  about  through 
careful  refinement  of  the  growth  process  to  the  point  where  pure  material 
can  be  achieved.  Then,  by  examining  the  low  temperature  (77K)  Hall 
mobility,  the  effects  of  alloy  scattering  can  be  separated  from  those 
due  to  ionized  impurity  scattering,  and  are  in  good  agreement  with  our 


previous  models. 
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Durlng  the  past  year,  one  phase  of  our  research  has  been  to 

investigate  other  quaternary  111-V  systems.  In  fact,  a complete 

survey  of  the  fifteen  major  quaternary  systems  of  the  form  y^y 

and  A^B^C^D  has  been  made.  The  systems  are  composed  of  the  compounds 

from  the  group  HI  atoms  Al,  Ga,  and  In  and  the  group  V atoms  P,  As,  and 

Sb.  Some  of  these  systems,  to  be  discussed  later,  also  have  desirable 

material  properties  for  electron  devices.  However,  we  have  taken  a 

wait-and-see  approach  before  submitting  this  data  for  publication  in 

the  open  literature,  since  we  feel  that  it  is  imperative  for  the 

calculations  on  Ga,  In  P,  As  to  be  verified  experimentally  to  some 
1-x  X 1-y  y 

extent  before  suggesting  other  now  quaternary  material  systems  for 
development.  Very  recently,  experimental  work  on  the  low  field  mobility 
of  the  ternary  Ga^^  ^In^As  lattice  matched  to  InP  has  been  reported  as  verifica- 
tion of  our  previous  calculations  for  this  material  [4].  It  should  be 
renujrked  tliat  much  progress  in  the  materials  growth  technology  has 
led  to  these  impressive  results. 

It  has  become  apparent  during  the  last  year  that  a significant 
emerging  research  area  is  that  of  submicron  device  physics  and  technology. 
During  the  past  grant  period  our  research  has  been  substantially  shifted 
to  emphasize  this  area  and  to  begin  a major  effort  in  this  field. 

We  would  like  to  concentrate  on  the  physical  processes  which  effect  the 
operation  of  small  devices  and  to  begin  to  use  the  Monte  Carlo  method 
for  this  purpose.  We  plan  to  question  based  on  first  principles  the 
physics  of  electron  transport  upon  which  the  Monte  Carlo  method  is 
based,  and  to  investigate  the  fundamental  limits  of  the  validity  of 


the  Monte  Carlo  method  as  a technique  for  studying  transport  phenomena 


in  general. 


For  example,  one  should  probably  begin  by  examining  a basic 
premise  of  most  traditional  treatments  of  transport  theory,  including 
the  Monte  Carlo  method,  which  separates  the  transport  process  into 
separate  and  non-interacting  processes  of  drift  in  the  field  and 
scattering  through  a lattice  or  defect  Interaction.  It  has  been 
pointed  out  that  under  the  conditions  of  high  electric  fields  and 
strong  scattering  rates  that  such  an  assumption  is  questionable  15], 
and  these  conditions  could  be  met  in  small  devices.  Our  plans  extend 
to  investigations  of  other  areas  including  statistical  concepts,  such  as 
fluctuation  phenomena,  modifications  of  collision  interactions  and 
interaction  time  field  effects,  velocity  overshoot  effects,  two-dimensional 
device  effects,  and  the  general  examination  of  device  concepts. 

There  is  one  final  point  to  be  discussed  here.  The  Monte  Carlo 
method  is  generally  considered  as  a single-particle  description  of 
transport  phenomena.  For  this  reason,  it  has  been  stated  that  the 
Monte  Carlo  method  will  have  difficulties  in  treating  small  device 
effects  because  electron-electron  interactions  cannot  be  conveniently 
included,  and  these  will  become  important  due  to  scaling  rules  which 
generally  require  higher  doping  levels  as  device  dimensions  decrease. 

Now,  besides  the  Monte  Carlo  method,  there  are  two  other  methods 
proposed  for  application  to  and  study  of  small  device  problems.  These 
methods  are  the  dr  if ted-Maxwell Ian  approximation  [6]  and  the  iterative 
integral  approach  [7].  Both  of  these  techniques  also  have  limitations 
at  high  doping  primarily  because  degenerate  statistics  are  difficult 
to  include  and  are  of  necessity  fundamentally  neglected  in  the 
formulations.  Also,  it  seems  to  us  to  be  premature  at  this  time  to 
argue  in  favor  of  the  dr  if ted-Maxwell ian  approximation  over  the  other  two 
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methods  as  being  the  high  electron  density  limit  to  the  distribution 
function  [8],  and  at  the  same  time  argue  in  favor  of  a reevaluation 
of  physical  concepts  as  applied  to  small  electron  devices.  For  example, 
the  drifted-Maxwelllan  may  be  a valid  velocity  distribution  function  at 
high  doping,  but  it  is  not  a valid  momentum  distribution  function  due  to 
non-parabolic  band  effects.  Also,  the  question  of  the  use  of  the  con- 
cept of  a distribution  function  must  be  raised  when  one  beings  to  consider  i 

small  devices  where  non-steady-state  processes,  such  as  velocity  over-  | 

shoot  and  time-dependent  diffusion  processes,  become  important.  The 
Monte  Carlo  method  does  not  necessarily  involve  the  concept  of  the 

distribution  function  except  in  a self-consistent  methodology,  and  , 

should  rightly  be  considered  as  the  major  fundamental  approach  to  modeling  I 

small  device  physics.  Our  approach  has  been  to  begin  to  apply  the  I 

Monte  Carlo  method  to  small  devices  by  making  major  modifications  to  the  I 

computational  techniques  which  are  oriented  toward  treating  such  small 
devices.  These  modifications  include  the  treatment  of  position-dependent 
electric  field  distributions  and  velocity-overshoot  and  diffusion-overshoot 
effects  in  general.  The  next  step,  which  has  been  accomplished  this  year,  | ^ 

is  to  also  include  degenerate  statistics  in  the  scattering  processes.  We 

believe  we  have  made  a major  contribution  in  this  regard,  and  papers  are  , 

being  prepared  for  publication.  j 

Thus,  a major  goal  of  this  research  is  to  examine  the  application  of 
the  Monte  Carlo  method  to  the  study  of  small  electron  devices.  It  is  strongly 

I 

believed  that  of  the  possible  techniques  to  be  appropriate  for  such  studies 

that  the  Monte  Carlo  method  will  again,  as  it  has  in  the  past,  prove  to  be  i 

the  more  complete  and  accurate  computational  method  to  study  transport  and 

hot  electron  phenomena  in  the  next  generation  of  electron  devices. 

I ^ 

> . 
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2.0  RESEARCH  RESULTS 

2.1  Review  of  Progress 

This  section  contains  a brief  discussion  of  work  performed  during 
the  past  year.  Also  Included  are  a list  of  publications  and  presentations, 
and  a summary  of  papers  in  progress  for  publication.  Finally,  the 
appendices  include  manuscripts  published  during  this  year. 

1)  The  calculation  of  iso-bandgap  and  iso-lattice  constant  contours 

for  Ill^-V^,  and  quaternary  alloys  has  been 

completed  and  published  during  1978  [9,10).  We  have  received 

several  requests  for  reprints,  and  a request  from  Rockwell 

International  for  additional  calculations  in  the  Al,  Ga  As,  Sb 

1-x  X 1-y  y 

quaternary. 

2)  A systematic  examination  of  the  transport  properties  of  fifteen 
quaternary  alloys  (nine  III2-V2  alloys  and  three  each  of  the 

and  alloys)  has  been  carried  out.  The  group  III 

elements  considered  are  Al,  Ga,  and  In  while  the  group  V elements 

are  P,  As,  and  Sb . Some  systems  have  properties  which  look  very 

promising,  especially  for  transferred  electron  devices.  None  of 

the  materials  seem  to  be  superior  to  Ga,  In  P,  As  for  FET's, 

' 1-x  X 1-y  y 

but  there  are  several  systems  possessing  high  peak  velocities, 

high  saturated  velocities,  and  low  threshold  fields  which  makes 

these  attractive  for  TED  devices.  For  example,  specific  systems 

are  Ga,  In  P,  Sb  , Ga,  In  As,  Sb  and 
1-x  X 1-y  y 1-x  X 1-y  y 


1 

4 


We  plan  to  submit  this  work  for  pvibl  Icat  ion,  but  would  like  to 

wait  for  more  advanced  materials  growth  technology  and  the 

measurement  of  material  arameters  for  Ga  In  P As  . 

1-x  X 1-y  y 

3)  Tiie  Monte  Carlo  program  has  been  modified  to  account  for  the  effects 

of  degeneracy  in  heavi ly-doped  materials.  At  low  electric  fields, 

wl\en  t tie  electron  ti'tnper.i  t ure  and  l.ittice  temperature  arc  equal, 

tlie  principle  modification  is  the  incorporation  in  the  scattering 

rates  of  an  occupancy  factor  l-t(K'),  wliere  E'  is  the  electron 

energy  .ifter  sc.atiering  and  t(E')  is  the  Eerml-Dirac  energy 

distribution  function  with  parameters  E,  and  T . At  low  fields, 

1 e 

these  parameters  can  be  determined  from  the  free  carrier  concentration 

3 

n and  the  average  carrier  energy  — kT^^.  At  high  electric  fields  the 
hot  electron  distribution  function,  .ind  thus  the  occupancy  factor, 
is  unknown,  but  the  Feimi-hirac  distribution  function  is  assumed 
to  be  a reasonable  first  approximation.  However,  in  this  case 
the  Fermi  energy  and  the  electron  lemiu’ralure  cannot  be  determined 
a priori,  since  the  electron  energy  is  nnknow'n.  The  approach  which 
we  have  taken  is  to  perform  the  Monte  Carlo  calculation  Iteratively, 
calculating  the  average  energy  in  each  iteration  and  using  this 

calculated  value  along  with  the  t'lectron  concent  r.at  ion  to  up-date 

t 

E-,  T , and  the  occupancy  factor  l-f(E).  The  Iterative  process 
i e 

l.s  continued  until  self-consistent  values  of  E^,  T , and  the 

t e 

.iverage  energy  are  obtained.  Those  values  are  then  used  in  the 
transport  calculations. 
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St'veral  other  modi  f Icat  tons  have  also  been  completed  in 
connection  with  the  simulation  of  heavl ly-dc ped  materials.  For 
example,  the  Inclusion  of  a transverse  magnetic  field  allows  the 
calculation  of  Hall  mobility,  and  the  required  modification  of 
Einstein's  relation  to  account  for  degeneracy  allows  the  accurate 
estimation  of  drift  velocity  at  low  applied  fields.  Some  of  the 
results  of  this  work  have  been  presented  at  the  Hot  Electron  Worksltop 
at  Cornell  Ibiiversity,  Aug.  17-18,  1978  and  at  the  7th  International 
Symposium  on  GaAs  and  Related  Materials  in  St.  Louis,  Sept.  2A-27, 

1978.  Two  jovirnal  articles  are  now  in  pt  ep.ira  t ion  . 

4)  The  Monte  Carlo  method  has  been  applied  to  the  determination  of 

the  spatial  distribution  function  and  its  temporal  evolution  for 

hot  electrons.  Among  other  things,  transients  in  carrier  velocity 

and  the  diffusion  process  itself  can  be  investigated  in  terms  of 

the  instantaneous  mean  m(t)  and  variance  s(t)  of  this  distribution. 

Velocity  overshoot  phenomena  such  as  obtained  by  other  techniques 

[11]  have  been  observed  and  can  be  stvidied  at  a more  basic  level  by 

this  approach.  We  have  also  observed  both  an  overshoot  and  an 

undershoot  in  the  diffusion  rate  or  diffusion  coefficient  D (here 

Id  2 

regarded  as  a function  of  time  through  tlie  definition  D(t)  “ T jY  ^ (t))- 
Transients  in  both  diffusion  coefficient  as  well  as  carrier  velocity 
persist  over  distances  traversed  by  carriers  which  .are  comparable  to 
important  dimensions  in  sub-micron  devices,  and  may  have  important 
implications  in  the  theory  and  design  of  such  devices.  Thus,  these 
transient  phenomena  warrant  further  study.  Preliminary  results  from 
this  work  were  presented  at  the  1978  Hot  Electron  Workshop  at  Cornell, 


■i 


and  a more  complete  discussion  will  be  contained  in  a paper  now 
in  preparation. 

5)  The  Monte  Carlo  program  is  currently  being  modified  to  allow  the 
addition  of  a sinusoidal  time-varying  electric  field.  The  basic 
approach  is  to  construct  a pseudo-ensemble  of  carrier  histories 
by  stacking  one-period  segments  from  a single  long-duration 
simulation.  This  modification  will  allow  Monte  Carlo  determination 
of  the  basic  frequency  response  characteristics  of  materials,  and 
should  provide  information  as  to  their  performance  at  mm-  and 
sub-mm  wave  frequencies.  This  work  is  in  the  initial  stages  and 

no  results  have  yet  been  obtained. 

6)  Several  other  areas  are  being  investigated  which  involve  major 
modif ica tions  in  the  Monte  Carlo  technique.  For  example,  the 
incorporation  of  boundaries  in  one-dimension  is  being  investigated 
in  order  to  allow  simulation  of  materials  with  finite  thickness. 

This  could  become  very  important  in  the  study  of  sub-micron  devices. 
Other  modifications,  which  will  be  discussed  in  the  next  section, 
are  the  incorporation  of  position-dependent  scattering  rates  and 
the  further  investigation  of  the  position-dependent  electric  field 
problem.  Also,  new  models  for  ionized-impurity  scattering  and 
electron-electron  scattering  which  examine  basic  statistical  factors 
in  the  fundamentals  of  these  problems  are  under  investigation. 

7)  Preliminary  work  on  the  transport  properties  of  II-IV-V2  chalcopyr ites 
has  been  completed.  Bulk  transport  properties  of  ZnSiAs^,  a material 
which  is  being  grown  in  our  laboratories,  have  been  studied  and 
certainly  look  very  interesting.  The  static  velocity-field  curve 
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r I'stMiil' 1 i‘S  tl\<it  nt  I III’  ill  iii.'iiiy  I t'.'ipi'i- 1 s . lli'wi'VtT,  tlu'if  .no  nui)i>r 
iiiu-oit  .'I  lilt  los  In  tlio  li.iiul  striH-tiiio  ot  tills  [lai  t lonlar  oi'iiii'oiind 
•ind  ot  t lu*  olia  1 1'opy  i 1 1 os  In  ^ionoi.il.  Wo  Irol  t h.it,  wo  i.lionld  lio 
oaiotnl  in  pi  osont  iiij^  sitoli  icsnlls  slnro  llu*y  .no  only  i.|ioonl  at  i vo 
at  pi  osi'iit  . Also,  tlioio  .no  m.i  )oi  iii.i  t ot  t .i  I s p.iowtli  pioMoms  in  tlio 
oli.i  1 ii'pv  I itos  wliioli  imisl  In'  solvoil  In'toio  t lioso  iiimponnds  will 
I'ooomo  iisoinl  oloot  ionic  iii.itoi  l.ils. 

K)  Work  on  t i .id  1 1 Iona  1 vo  1 oc  1 1 y ovoi  shoot  c.i  1 v' u 1 a t I ons  (11)  h.is  hoon 

complotod  lor  Oa,  In  1’,  As  . I'lio  losnils  snp^ost  that  it  liipli- 

1 - X X 1 y y 

mohilits'  ui.iloi  i.'il  with  a h.iiui-n.ip  ol  ahont  1 »'V  could  ho  p.i  own  on 
Ini'  snhslt.it  os,  thou  volocltv  overshoot  ol  tools  t.lionld  ho  ohsoiv.ihlo 
in  dovicos  wllli  I'h.innol  loiipths  holwoon  0.  “i  iiiii  - I . ()  |iin.  Tlio.  oxioiil 
ot  t ho  ovi-ishoot  is  .snl>sl  .1111  i.i  1 1 V pio.ilor  I h.in  Ih.il  tor  ('..lAs  i>i  liiT 
in  this  ipi.it  I'l  nary  . I'l  hoi  ipi.i  t oi  nai  i os  .ind  I oi  n.n  ioi.  at  o hoin^ 
ox. nil  1 noil  .It  this  I imo . 

‘1 ) A I Wi>-<l  i iiK'iis  i on.i  1 dovioo  model  li.is  hoon  liovolopod  .mil  is  now  lu'iii^; 

liirlhoi  I ol  i nod . Tho  oompnl.it  ion.il  piooodiiio  is  voiv  simil.n  to  Ih.il 
hoinp,  used  hv  H.  1..  llinl'in  at  linilod  Too  lino  1 op  I os  Kio-.o.i  j cli  I'onloi. 

Tho  mod  i I leal  ions  now  hoinp  1 ncoi  |'o  i .i  I od  involve  tlio  nso  ol  tlio 
ipuis  i - 1 i no.n  1 .M I i on  toolmiipio  which  h.ts  hoon  ntili/.od  N.  t'.  St  .it  o 

Ihilvotsitv  to  iiH'dol  i.i'l.n  I'olls  | 1 | . Wi-  hi’lli'vo  this  loolnilipio  will 

onh.iiioo  the  i.ii'.ih  1 1 i t I os  ol  tho  dovioo  piopi.nn.  I'lio  lonp  loim  p.oal 
is  to  piovlilo  .1  dovioo  piop.i.im  whiih  will  intoi.iot  with  tlio  Monte 
I'.nlo  piopi.nn  to  piovldo  a soil  oonstslont  '.olntioii  lot  I i. no. poll 
and  loimin.il  pi  opot  t los  ol  i.iih  mioiini  dov  i.  os  . 
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10)  Work  has  been  carried  out  to  experimentally  and  theoretically  study 

alloy  scattering  In  Ga.  In  P.  As  . The  experimental  work  has  been 

i ""X  X y 

supported  by  AFOSR  and  was  done  in  conjunction  with  the  Research 
Triangle  Institute.  The  results  were  presented  at  the  7th  International 
Symposium  on  GaAs  and  Related  Compounds  held  in  St.  Louis,  Sept.  2A-27, 
1978,  and  they  appear  to  suggest  that  alloy  scattering,  or  some  other 
form  of  defect  scattering  having  a temperature  dependence  similar 
to  alloy  scattering,  is  greatly  affecting  the  low-field  mobility 
in  this  quaternary.  One  such  mechanism  could  be  space-charge  or 
precipitate  scattering  [13].  This  point  is  made  in  light  of  the 
results  presented  by  Pearsall  [4]  at  the  same  meeting  on  Ga  ^.^In  53^® 
(lattice  matched  to  InP)  which  were  interpreted  to  suggest  that  alloy 
scattering  was  not  significant  (at  room  temperature)  in  this  ternary. 

It  is  very  probable  that  future  advances  in  growth  technology  for 

Ga^  ^'’x^l  y^^  will  lead  to  substantial  improvements  in  the  transport 

properties. 

11)  The  work  on  negative  resistance  in  the  central  1(000)  valley  of 
compound  semiconductors  has  been  completed  and  accepted  for 
publication  in  Solid-State  Electronics.  The  effect  is  real  a.ad 

our  previous  interpretation  of  it  stands  as  originally  reported.  The 
inclusion  of  the  p-state  mixing  phenomena  resulted  in  no  changes 

in  the  original  conclusions,  and  the  phenomenon  appears  to  have 
other  previously  unforseen  implications.  For  example,  velocity 
overshoot  in  materials  whose  peak  velocity  is  determined  by 


I 
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t'«'ntral  valloy  si-at  t or  lug  Is  illllfn-nt  from  that  of  matorfals 
such  as  CaAs  and  Ini’.  The  ov«*rshi>ot  soi-ms  to  persist  for  much 
longer  l iiiu's  In  such  materials.  Thus,  ft  ts  piissfhle  that 
central  va  1 1 ey -ilomlnat  ed  compv'und  .sem  Ici'iuluc  t vM  s could  he  v»‘iy 
tmpiutant  In  svih-mlcron  device  app  1 f cat  ions  . 
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CARRIER  COMPENSATION  AND  ALEOY  SCATTERING 

IN  Ga,  In  P,  As  GROWN  BY  LIQUID  PHASE  EPITAXY 
1-x  X l-y  y 

M.  A.  Littlejohn,  R.  A.  Sadler,  T.  U.  Glisson,  and  J.  R.  Hauser 
Electrical  Engineering  Department 
North  Carolina  State  University 
Raleigh,  North  Carolina 


ABSTRACT 


Carrier  mobilities  achieved  experimentally  in  Ga,  In  P,  As 

^ ^ 1-x  X 1-y  y 

alloys  are  considerably  less  than  those  predicted  by  Monte  Carlo  transport 
calculations.  In  this  paper,  a comparison  between  the  temperature 
dependence  of  electron  Hall  mobility  measured  experimentally  on  LPE- 
grown  samples  and  calculated  from  a Monte  Carlo  method  is  presented. 

The  experimental  data  is  obtained  from  samples  which  are  n-type  and 
are  grown  on  semi -insulating  InP  substrates.  The  theoretical  calculations 
are  based  on  a diffusion  coefficient  estimator  for  calculation  of  Hall 
mobility.  An  analysis  of  the  data  indicates  that  both  alloy  scattering 
and  carrier  compensation  are  necessary  to  explain  the  observed  temperature 
dependence  of  Hall  mobility. 


1.  Introduction 


Tlie  III-V  quaternary  alloy  Ga,  In  P,  As  has  received 
considerable  attention  recently  due  to  many  desirable  electronic  and 
material  properties,  and  this  compound  has  been  used  in  a wide  variety 
of  device-related  studies  (Clawson  1978).  In  spite  of  the  wide 
variety  of  developed  applications,  the  desirable  transport  properties 
proposed  for  this  material  based  on  calculations  using  the  Monte  Carlo 
method  (Littlejohn  e^  aJL  1977a)  have  not  been  fully  realized  (Moiioc  et  al . 
1978,  Houston  1978).  In  particular,  the  reported  Hall  mobilities 

have  been  considerably  lower  than  predicted,  and  room  temperature 
electron  concentrations  below  about  IxlO^^cm  ^ have  been  difficult  to 
achieve.  Much  of  the  work  has  been  done  using  liquid-phase  epitaxy 
(LPE).  However,  the  electron  concentrations  are  much  larger  for  the 
quaternary  than  those  usually  achieved  for  LPE  growth  of  the  binary 
compounds  GaAs  and  InP,  even  though  the  starting  constituent  source 
materials  used  to  grow  the  quaternary  are  of  the  s.une  purity  as  those 
used  for  growth  of  the  binaries. 

In  this  paper  a comparison  between  experimental  Hall  mobilities  of 
LPE-grown  quaternary  layers  on  InP  substrates  and  those  calculated  from 
the  Monte  Carlo  method  is  presented.  Based  on  this  comparison  it  is 
concluded  that  both  carrier  comju'nsat i on  as  well  as  alloy  scattering 
are  contributing  to  the  low-field  Hall  mobility  and  to  its  temperature 
dependence . 
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2.  Experimental  Aspects 

The  samples  studied  in  this  paper  were  grown  by  liquid  phase 
epitaxy  at  660°C  in  a standard  horizontal  slider  system  (Phatak  and 
Harrison  1978).  The  as-grown  layers  were  undoped,  and  the  substrate 
material  was  semi-insulating  Fe-doped  InP  in  a (lll)B  orientation. 

Typical  layer  thicknesses  were  around  10  pm. 

The  Hall  measurements  were  made  in  a variable  temperature  cryostat 
between  77K  and  300K.  A standard  Van  der  Pauw  sample  geometry  was  used, 
and  contacts  to  the  samples  were  made  with  an  In-Sn  alloy  (Houston  £t_  al_. 
1978). 

3.  Monte  Carlo  Calculations 

The  calculations  presented  in  this  paper  were  obtained  by  the 

Monte  Carlo  method,  using  a program  described  previously  (Littlejohn 

et  al . 1977a).  Of  particular  importance  is  the  method  used  to  obtain 

the  quaternary  material  parameters  for  Ga^  x^'^x^l  done 

by  an  interpolation  procedure  which  is  based  on  an  assumed  knowledge  of 

the  constituent  binary  material  parameters  of,  in  this  case,  GaAs,  InAs, 

InP,  and  GaP.  IVhile  much  needs  to  be  done  to  further  validate  this 

interpolation  procedure,  it  has  proved  useful  in  calculating  energy  band 

gap  and  lattice  constant  (Glisson  e;^  al_.  1978,  Moon  e^  a^.  1974)  as  well 

as  electron  effective  mass  (Restorff  £t  al_.  1978).  In  these  cases  it  has 

led  to  good  agreement  with  experimental  data.  Figure  1 shows  the  results 

of  a calculation  for  the  T-valley  electron  effective  mass  in  Ga,  In  P,  As 

1-x  X 1-y 

as  a function  of  As  composition  for  layers  lattice-matched  to  InP.  The 
three  different  curves  represent  calculations  based  on  different  values 


I 
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of  binary  material  parameters.  However,  for  y>0.5  the  effective 

I mass  does  not  strongly  depend  on  the  uncertainty  in  the  InP  effective 

I mass  (y=0) (Maloney  and  Frey  1977). 

In  general,  the  Monte  Carlo  method  is  difficult  to  apply  with 
good  statistical  accuracy  at  low  fields  and  low  temperatures.  However, 
this  problem  can  be  overcome  to  a large  extent  by  using  statistical 
estimators  for  the  diffusion  coefficients  (Canal  i ^ al_.  1975). 

In  addition,  the  present  program  has  been  further  modified  to  include 
the  effects  of  heavy  doping  as  well  as  to  use  crossed  electric  and 
magnetic  fields  for  the  computation  of  Hall  mobility,  since  it  is 
Hall  mobility  and  not  drift  mobility  that  is  usually  measured.  The 
program  details  will  be  presented  in  a future  publication.  Heavy  doping 
effects  can  be  important  in  Ga,  In  P,  As  alloys  having  small  effective 
masses  such  as  shown  in  Figure  1.  For  GaAs,  with  m*=0.063  m^,  the 
conduction  band  effective  density  of  states  is  about  4xl0^^cm  ^ at  300K. 
From  the  values  in  Figure  1 a number  of  1.65xl0^^cm  ^ for  the  density 
of  states  of  a quaternary  with  y=0.7  is  calculated.  As  the  temperature 
is  decreased  to  77K,  the  conduction  band  effective  density  of  states  is 
reduced  to  SxlO^^cm  ^ for  GaAs  and  to  Z.lxlO^^cm  ^ for  the  quaternary. 

Thus  one  could  expect  degenerate  statistics  to  be  important  in  determining 
the  quaternary  transport  properites  for  presently  achievable  doping 
1 evcls. 

The  Nionte  Carlo  method  has  the  additional  feature  that  very  detailed 
transport  and  scattering  factors  can  he  included  in  the  calculations. 
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These  factors  include  such  band  structure  effects  as  p-state  mixing 
and  band  non-parabol icity  (Fawcett  £l_-  1970).  Also,  the  Monte  Carlo 
method  takes  into  account  thermal  stimulation  at  temperatures  above  300K 
of  carriers  into  higher-lying  conduction  band  minima  (Rode  1975) . If 
care  is  taken  in  assuring  statistically  accurate  estimates  in  the  Monte 
Carlo  method,  then  the  technique  should  be  accurate  and  very  complete, 
and  as  convenient  as  techniques  based  on  classical  transport  analysis 
and  the  use  of  other  numerical  methods  (Rode  1975). 

4 . Discuss  ion 

The  use  of  the  Monte  Carlo  method  for  calculating  Hall  mobility  was 
first  examined  by  considering  data  for  GaAs  since  it  was  felt  that 
values  for  all  important  transport  coefficients  and  material  parameters 
were  reasonably  well  known  (Rode  1975,  Littlejohn  1977b, 

Kratzer  and  Frey  1978).  Figure  2 shows  Hall  mobility  as  a function  of 
free  electron  concentation  for  GaAs  at  300K  for  a compensation  ratio  of 
l(i.e.  N +N  /n  = 1 , where  N*  = ionized  donor  density,  N = ionized 
acceptor  density  and  n = free  electron  density).  ShowTi  for  comparison  in 
this  figure  is  a range  of  experimental  data,  taken  for  convenience  from  the 
last  two  proceedings  of  this  conference  and  the  standard  data  from  Sze,  19(i9. 
Care  was  taken  to  record  only  data  reported  as  Hall  mobility  data.  Also 
shown  are  the  results  of  a calculation  for  Hall  mobility  based  on  an 
iterative  integral  method  (Rode  1975).  Tliere  are  some  discrepancies 
between  these  two  approaches,  as  can  be  seen  in  this  figure.  For  low 
carrier  concentrations  the  discrepancies  (about  13i>)  can  be  attributed 


l; 
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to  slightly  different  values  of  bulk  material  parameters  used  in  the  two 

calculation  methods.  The  material  parameters  used  here  in  the  Monte  Carlo 

calculation  have  been  chosen  based  on  both  low-field  mobility  values  and  on 

high  field  transport  properties  (Littlejohn  1977b,  Kratzer  and 

Frey  1978).  IkDwever,  the  Monte  Carlo  method  will  yield  Hall  mobilities 

14-3  2 

at  low  carrier  concentrations  (<10*  cm  ) well  above  9000  cm  /volt  sec 
if  slight  adjustments  in  bulk  parameters  are  made,  and  certainly  the 
sensitivity  of  the  Hall  mobility  in  pure  GaAs  to  slight  variations  in  material 
parameters  is  well-recognized  (Rode  1975).  At  high  carrier  concentrations 
(>10^^cm  ^)  the  deviations  between  these  two  computational  methods  are  much 
more  subtle.  They  are  due  to  the  way  in  which  the  ionized  impurity 
scattering  rate  is  corrected  for  degeneracy  in  the  Monte  Carlo  method  and 
this  will  be  a subject  for  a future  paper.  However,  in  spite  of  the 
differences  the  authors  believe  that  the  Monte  Carlo  method  is  an  accurate 
technique  for  calculating  Hall  mobility  which  gives  agreement  with 
experimental  data  over  a wide  range  of  carrier  concentrations. 

Figure  3 shows  experimental  Hall  mobility  data  versus  temperature 
for  pure  GaAs,  along  with  the  calculated  temperature  dependence  of 
mobility  based  on  both  the  iterative  integral  and  the  Monte  Carlo  method. 

Here  the  Monte  Carlo  method  gives  an  excellent  fit  over  the  temperature 
range  150K-600K,  while  Rode  has  pointed  out  that  the  iterative  integral 
lies  about  10%  above  the  experimental  data  in  this  range  (Rode  1975). 

IVhi'e  we  have  not  done  so  in  this  study,  it  should  be  pointed  out  that 
the  Monte  Carlo  method  can  be  extended  to  temperatures  below  lOK 
(Canali  1975).  Also,  other  GaAs  material  has  been  successfully 

studied  by  the  Monte  Carlo  method,  and  the  general  use  of  the  method 
appears  to  be  well-justified. 


Shown  in  Figure  2 is  a calculation  of  the  Hall  mobility  versus 

electron  concentration  at  300K  for  Ga  „^In  _,P  .As  , based  on  the 

.27  .73  .4  .6 

Monte  Carlo  method  and  the  material  parameter  interpolation  procedure 

described  previously  (Littlejohn  al_.  1977a).  Presently,  to  the 

authors'  knowledge,  the  largest  value  of  room  temperature  Hall  mobility 

2 

for  any  quaternary  is  about  6000  cm  /volt  sec  with  an  electron  . 

concentration  of  about  7xl0^^cm  ^ (Houston  e^  al . 1978).  Figure  4 

shows  a summary  of  typical  Hall  mobility  versus  temperature  for  the 

quaternary  (and  the  InP  lattice-matched  ternary  Ga  j^^s)  taken  in 

our  laboratory  and  reported  in  the  literature  (Houston  ^ al_.  1975, 

Takeda  al_.  1976).  In  general,  the  room  temperature  Hall  mobilities 

2 

seem  to  lie  in  the  range  of  3000-4000  cm  /volt  sec  with  electron 

16  “3 

concentrations  of  (2-4)xl0  cm  . Interestingly  enough,  the  highest 
mobilities  and  lowest  electron  concentrations  consistently  seem  to  occur 
near  the  ternary  boundaries  for  the  quaternary,  i.e.  for  either  relatively 
low  As  or  P compositions.  This  is  suggestive  of  the  influence  of  alloy 
scattering  (Littlejohn  al_-  1978). 

To  investigate  the  general  problem  of  carrier  compensation  and  alloy 
scattering  in  the  quaternary  Ga^  x^^x^l  y^^y  Motite  Carlo  method  was 
used  to  study  several  typical  samples.  Figure  5 illustrates  the  situation 
when  an  attempt  was  made  to  fit  the  data  using  either  carrier  compensation 
with  no  alloy  scattering  or  alloy  scattering  with  no  carrier  compensation. 
First,  the  alloy  scattering  potential  (Littlejohn  e^  al . 1978)  was 
varied  while  keeping  the  compensation  ratio  equal  to  1 (i.e.  N =N  =n) 
until  agreement  was  obtained  between  the  measured  and  calculated  Hall 
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mobility  at  room  temperature.  Then  the  temperature  dependence  of  Hall 
mobility  was  calculated,  resulting  in  curve  a)  of  Figure  5 which  is  in 
very  poor  agreement  with  the  data.  Curve  b)  results  by  removing  alloy 
scattering  from  the  calculation.  First,  the  ionized  impurity  density  was 
varied  until  agreement  was  achieved  with  the  measured  Hall  mobility  at 
room  temperature,  while  keeping  the  free  carrier  density  equal  to  its 
measured  room  temperature  value.  Curve  b)  is  the  resulting  temperature 
variation  of  Hall  mobility  with  only  carrier  compensation  present. 

Again  the  agreement  with  experimental  values  at  other  temperatures  is 
very  poor,  and  it  was  concluded  that  neither  alloy  scattering  nor  carrier 
compensation  could  account  for  the  measured  temperature  dependence  of 
Hall  mobility. 

Next,  an  attempt  was  made  to  find  a consistent  set  of  values  for 
both  the  alloy  scattering  potential  and  the  ionized  impurity  density  which 
would  together  give  a good  fit  to  the  temperature  dependence  of  Hall 
mobility.  To  do  this  two  experimental  points  were  chosen,  one  at  300K 
and  one  at  about  150K.  A detailed  parameter  study  was  made  by  varying 
the  alloy  scattering  potential  and  ionized  impurity  density  until  a set 
of  values  for  these  two  parameters  which  agreed  with  the  measured  values 
of  both  Hall  mobility  and  electron  concentration  at  the  two  temperature 
extremes.  Then  the  complete  Hall  mobility  versus  temperature  curve  was 
calculated  using  these  values.  The  results  are  shown  in  Figure  6.  The 
Monte  Carlo  estimates  lie  within  the  dashed  curves  (several  estimates 
were  made  at  each  temperature)  while  the  bars  indicate  a 5%  variation 
in  the  experimental  data  at  the  two  temperature  extremes.  It  should  be 


pointed  out  that  the  fits  are  facilitated  by  the  fact  that  the  electron 
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lomciit  rat  ion  lUn-s  not  vary  s i>;in  f i rant  1 y for  tlu'so  niiat  i*rnary  sampli's 
bctwcon  100K-;U)0K.  The  compensation  ratio  of  3.9  is  very  close  to  a 
sample  with  a very  similar  temperature  ilepeiulence  of  mobility  studied 
hy  Houston  1978. 

I’lio  calculations  are  more  sensitive  to  the  ionised  impurity  density 
tlian  t))ey  are  to  the  alloy  scattering  potential.  Thus,  it  is  not 
feasible  to  ilraw  many  conclusions  from  the  .actual  value  ol  the  alloy 
scattering  parameter  shown  in  figure  li.  The  value  is  larger  tlian  the 
prciposed  theoretical  values  for  three  ilifferent  moilels  (l.iltlejohn 
et  aj_.  1978)  by  as  much  as  iO*!,  to  lOOI.  Tliere  is  another  scattering 
mech.'inism  which  has  the  same  temperature  ilei'eiulence  of  mobility  as  alloy 
scattering,  which  is  often  called  space  charge  scattering  (Weisberg  19<i2). 

If  this  mechanism  were  also  present  in  these  samples  it  would  be  dilticult 
to  detect  from  il.ill  measurements  alone,  but  could  ellect  the  determination 
of  alloy  scattering  luitential  by  the  teciinique  used  here.  This  added 
effect  could  be  the  result  of  prec  i p i t .it  es , clustering,  or  other  types 
of  second  ph.ise  defects. 

Other  samples  with  similar  Hall  mobility  temiierature  dependences 
have  been  studied  with  ver>’  similar  ri'sults  as  given  in  figure  0.  In 
giMier.il  for  these  qu.it  ern.i  r i es  the  weaker  the  t emper.it  ure  depeiulcnce  ot 
Hall  mobility  between  .ibout  I OOk-.lOOK,  the  more  important  alloy  scattering 
lu'ciniies . 

It  is  interesting  to  considi'f  those  s.imi'les  with  a stronger  teiiq'erature 

dependence  in  this  t eni]>erat  ure  range,  figure  ’ shows  the  results  (or  a 

C.a  pin  ^,/\s  samjile,  which  is  a ternary  in  this  qua  t ern.i  ry  system  having 

•> 

.1  Il.ill  mobility  at  room  t <'mjier.il  ure  l.irger  than  8(100  cm  /volt  sec  (Morkoc 
1978,  T.ikeda  1978).  for  this  case,  the  temin'rature  ilepeiulence  of  Il.ill 


w 
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inobility  could  be  fit  by  the  Monte  Carlo  method,  but  with  an  alloy 
scattering  parameter  essentially  the  same  as  that  predicted  by  the 
theoretical  model.  Tlie  influence  of  alloy  scattering  in  this  ternary 
is  less  imiH>rtnnt  in  determining  the  Hall  mobility  than  it  is  in  the 
quaternary.  Tins  fact  is  somewhat  surprising  based  on  previous  results 
for  the  scattering  potentials  in  both  ternary  and  qiuUernary  alloys 
(Littlejohn  ^ 1978). 

.5.  Conclusions 

Alloy  scattering  appears  to  play  an  imjH'rtant  role  in  the  low-field 

transport  properties  of  C>a.  In  l\  As  . This  conclusion  is  baseil  on  the 
' ‘ * 1-x  X 1-y  y 

tlata  and  calculations  jiresented  in  this  iKqier,  as  well  as  on  other 
published  Hall  effect  data.  However,  the  accuracy  of  the  calculations 
are  preilicated  on  the  material  parameter  int  erj'olat  ion  procedure,  which 
is  presently  required  for  any  transport  calculations  in  this  quaternar\. 
Historically,  the  Monte  (^irlo  method  has  become  more  useful  in  studying 
transjxirt  jiropert  ies  in  compouiul  semiconductors  as  accurate  material 
parameters  become  available.  I'he  ileterminat  ion  of  the  basic  material 
parameters  should  be  an  imixirtant  part  of  programs  which  attempt  to 
obtain  high-j>urity  Caij  x'*'x''l  y^^y  H"'*' *'*'”'* 'V  alloys. 
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Llst  of  Figures 


Figure  1.  Effective  mass  ratio  as  a function  of  As  composition  In 

Ga,  In  P,  As  . Curve  (a)  uses  binary  effective  masses 
1-x  X l^y  y 

from  Restorff  1978.  Curves  (b)  and  (c)  use  values  from 
Littlejohn  1977a,  with  two  different  reported  values  for  InP. 


Figure  2.  Hall  mobility  versus  free  carrier  concentration  for  a 

compensation  ratio  of  1 . The  legend  Is:  ———Monte  Carlo 
calculations,  Iterative  integral  calculations  (Rode  1975), 

Data  from  Sze  1969,  © - Other  experimental  data. 


Figure  3.  Hall  mobility  versus  temperature.  The  legend  is: 

- Monte  Carlo  calculations,  and  © - Iterative 

integral  calculations  and  experimental  data  from  Rode  1975. 


Figure  A.  Typical  Hall  mobility  data  versus  temperature  for  Gaj^_^In^Pj^  As  . 
/ The  curves  show  data  from  Houston  et  1978, 

Takedl  e£  al.  1976,  and  from  this  study. 

Figure  5.  Hall  mobility  versus  temperature  for  Ga  ^In  .46^®. 54' 

The  dotted  points  are  experimental  data  while  curves  a)  and 
b)  represent  Monte  Carlo  calculations  with  no  carrier 
compensation  and  no  alloy  scattering,  respectively. 

Figure  6.  Experimental  data  and  Monte  Carlo  calculations  of  Hall  mobility 
versus  temperature  for  Ga  ^In  ^^As  showing  fitted 

material  parameters. 


Figure  7.  Experimental  data  and  Monte  Carlo  calculations  of  Hall  mobility 
versus  temperature  for  Ga  53^®  showing  fitted  material 

parameters. 
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ENERGY  BANDGAF  AND  LATTICE  CONSTANT 
CONTOURS  or  in-v  quaternary  alloys* 


T.  H.  Glisson,  J.  R.  Kauser,  M.  A.  Littlejohn,  C.  R.  Williams 

North  Carolina  State  University 
Electrical  Engineering  Department 

Raleigh,  North  Carolina  27607 


(Receivixl  May  Ih.  1*^77) 


Energy  band  gap  and  lattice  constant  contovirs  are  presented 
foi'  the  nine  i^uaternar\’  .tlloys  lormed  trem  Al,  ^a , In  and 
P,  As,  Sb.  The  quaternary  handgaps  were  obtained  using  an 
interpolat  ion  forntila  proposed  by  'loon  7he  quater- 

nary lattice  constants  were  obtained  by  use  of  a linear 
interpolation  technique  using  the  bin. try  lattice  constants 
as  boundary  values. 

Key  words:  qv:atein.try  alloys,  bandgap,  lattice  constant. 


I n t ivkI uct  ion 

There  has  been  considerable  interest  in  the  quater- 
n.jry  III-V  semiconductor  materials  lor  many  .ippl i cat  ions  , 
such  as  e lect ro -op t i cs  and  microwave  devices  [l-6j.  One 
significant  reason  for  this  interest  is  the  ability  to 
synthesize  materials  with  a fixed  lattice  constant  and  a 
variable  range  of  energy  handgaps  or,  cor.versely,  to 
synthesize  materials  with  a constant  energy  bandgap  and 
a variable  lattice  constant  [7]. 

■•This  work  was  supported  by  a research  grant  fi'om  the 
Office  of  Naval  Research,  Washington,  D.C. 
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In  order  to  predict  the  energy  bandgap  or  lattice 
constant  for  a quaternary  material  a two-step  sequence  is 
involved.  Initially,  it  is  necessary  to  have  available 
experimental  (or  theoretical)  values  of  these  material 
parameters  for  the  four  possible  binary  III-V  constituents 
of  the  quaternary  compound.  These  values  are  readily 
available  for  most  III-V  materials  in  the  literature  [8]. 
There  are  also  four  possible  ternary  (pseudobinary)  III-V 
systems  in  a given  quaternary  compound.  In  the  first  step, 
the  energy  bandgap  and  lattice  constant  as  a function  of 
alloy  composition  for  each  ternary  system  is  computed 
from  the  known  binary  material  parameters.  This  computa- 
tion is  on  a good  experimental  and  theoretical  basis 
[9-13]  involving  the  use  of  Vegard's  law  and  the  concept 
of  the  "bowing  parameter"  for  the  energy  bandgap  [1**]. 

The  final  step  involves  the  use  of  the  compositional 
dependence  of  the  energy  bandgap  and  lattice  constant  for 
the  four  ternary  combinations  to  interpolate  to  the  desired 
quaternary  energy  bandgap  and  lattice  constant.  This  final 
interpolation  step  is  somewhat  empirical  in  nature.  There 
appear  to  be  several  such  interpolation  techniques  used  in 
the  literature  [6,7,15]  and  while  these  are  similar,  they 
also  have  non-trivial  differences.  This  paper  will  discuss 
this  critical  interpolation  procedure  in  view  of  the  small 
amount  of  experimental  quaternary  data,  and  will  present 
calculations  for  energy  bandgaps  and  lattice  constants 
for  the  nine  quaternary  alloys  formed  from  Al,  Ga , In  and 
P,  As,  Sb  using  one  of  these  interpolation  techniques. 


Notation 


There  appears  to  be  no  generally-accepted  scheme  for 
symbolically  describing  the  III-V  quaternary  alloys.  In 
this  paper  the  following  convention  has  been  adopted. 

Within  both  the  group  III  and  group  V pairing,  the  first 
position  is  occupied  by  the  element  with  lowest  atomic 
number.  The  composition  variables  are  associated  with  the 
four  elements  in  the  order  1-x,  x,  1-y,  y.  for  example,  a 
quaternary  alloy  will  be  denoted  as  A^  x®x^l  ^ 

and  B are  group  III  elements  with  A having  lower  atomic 
number  than  B.  Likewise,  C and  D -n'e  group  V elemiuits  with 
C having  lower  atomic  number  than  D.  Following  previous 
authors  [6,7],  a quaternary  alloy  parameter  (e.g.,  bandgap 
or  lattice  constant)  is  described  by  a surface  Q(x,y)  over 
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the  x,y  composition  plane  (0<  x < 1,  0 < y < 1).  At  the 
corners  (x  and  y equal  to  zero  and/or  one)  the  values  of 
the  parameter  for  the  four  binary  elements  Q(0,0)  = , 

Q(1,0)  = B , Q(l,l)  = B^  and  Q(0,1)  = are  obtained. 
Along  the  Boundaries  of  th«“  plane  the  parameter  for  the 
four  ternary  ele-ients  Q(x,0)  = 

Q(x,l)  = Tj^jCx)  and  0(0, y)  = obtained. 


Estimatioi.  of  Alloy  Pararioters 

Many  ternary  alloy  parameters  (e.g.  lattice  constant) 
are  obtained  by  linear  interpolation  from  those  of  the 
constituent  binary  compounds,  i.e., 

T.  .(x)  = xB.  t (l-x)B.  (1) 

For  some  parameters  the  theoretical  variation  with  compo- 
sition is  nonlinear  [14,16].  For  the  direct  bandgap  for 
example,  Thompson  and  Woolley  [9]  have  shown  that 


T..(x)  = xB . + (l-x)B.  - C. .x(l-x) 
13  3 1 13 


(2) 


where  C.  . is  the  bowing  parameter  for  the  ternary  alloy 
bandgap 

In  the  absence  of  definitive  theories  for  quaternary 
parameters,  estimates  of  a quaternary  alloy  parameter 
Q(x,y)  must  be  obtained  by  interpolation  from  the  four 
ternary  alloy  parameters  T...  Various  interpolation 
schemes  have  been  proposedi^^  Onton  and  Chicotka  [15] 
used  the  solution  of  Laplace's  equation  subject  to  the 
boundary  conditions  Q(x,0)  = etc.  In  Monte  Carlo 

transport  studies  [6]  the  present  authors  have  used  the 
interpolation  equation 


Q(x,y)= 


x(l-x)[(l-y)T^2(x)tyT^^(x)]+y(l-y)[(l-x)T^^(y)+xT^^(y)] 




x(l-x)+y( l-y ) 


14' 


23 


(3) 


This  int- . ,'>olat  ion  equation  reduces  to  the  ternary  param- 
eters on  the  quaternary  plane  boundaries  and  to  the  average 
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of  the  ternary  parameters  at  the  midpoint  (x=0.5,y=0.5) 
of  the  compositional  plane.  This  is  incidentally  the 
solution  to  Laplace's  equation  if  the  T. . are  linear  as 
in  Eq.  (1).  Moon  et  a^.  [7]  have  propoi^d  similar  schemes 
for  the  estimation  of  lattice  constant  and  bandgap.  For 
the  lattice  constant,  Moon  gives  (in  the  above  notation) 


For  the  lattice  constant,  the  T. . are  assumed  to  be 
linear  in  their  argumi>nts,  so  the  soiStion  to  Laplace's 
equation  is  identical  to  Eq,  (3),  which  in  turn  also 
reduces  to  Eq.  (4). 


For  the  bandgap.  Moon  e^  al^.  have  used  the  equation 


where  the  T. . are  determined  from  Eq . (1)  and 


An  important  difference  between  tliis  bandgap  estimate  and 
that  obtained  from  Eq.  (3)  is  in  tlie  iiunner  in  which  the 
ternary  bowing  parameters  C..  enter  the  calculation.  In 
Eq . (3),  the  ternary  bowing^t^f  fects  are  included  in  the 
T.  whereas  in  Moon's  nK!thod  the  bowing  enters  as  a 
separate  quaternary  bowing  pai'ametei'  term  given  in  Eq . (6). 
At  the  center  of  the  composition  pi. me,  Eq.  (3)  yields 
(for  the  quaternary  bowing)  1/16  + ^^3^ 

and  Eq.  (5)  gives  1/8  (C^^j  + ^-’3^’ 

twice  as  lai-ge.  Tims,  at  the  center,  Eq.  t3)  gives  the 
aver.age  of  the  bowing  contribut  ions  from  the  four  ternaries, 
whereas  Eq.  (5)  gives  the  average  of  the  bowing  conti'ibu- 
tions  from  the  two  sublattices.  This  is  the  principle 
difference  between  the  interpolat ion  schemes  of  Refs.  6 
and  7 . 


I 

I 

I 

I 

( 

I 

! 


-39- 


Encrgy  Bandgap  of  lll-V  Qiialemary  Alloyi 

Comparison  with  Experiment 
Onton  and  Chicotka  compared  their  interpolation 


approach 

(solution  to 

Laplace ' s 

equation)  with 

measured 

values  of 

the  quaternary  bandgap  in  Ga^  x^^x^l 

As  and 

found  the 

interpolated  values  to  be  within  +30 

mev  of  the 

measured 

values.  We 

have  repeated  this  comparison,  using 

their  measured  data. 

The  results  are  given  in 

Table  1. 

The  average  error  and 

standard 

deviation  of  the 

? error 

are  also 

shown,  and  it  is  seen 

that  Eqs.  (3)  and  (5)  give 

comparable  errors. 

Table  I. 

Comparison 

of  Eqs  (3) 

and  (5)  with  Experimental 

Data 

for  Ga, 

In  P.  As  . 

1-x 

x 1-y  y 

Composition 

Bandgap  (eV) 

X 

y 

Eq  (3) 

Eq  (5)  Measured  [15] 

0.004 

0.260 

2.357 

2.353 

2.222 

0.006 

0.340 

2.244 

2.239 

2.  124 

0.006 

0.420 

2.136 

2.131 

1.993 

0.009 

0.500 

2.027 

2.020 

1.881 

0.008 

0.590 

1.913 

1.907 

1.738 

0.010 

0.710 

1.763 

1.755 

1.592 

0.012 

0.830 

1.619 

1.609 

1.475 

O.OIV 

0.910 

1.521 

1.510 

1.351 

0.015 

0.970 

1.455 

1.447 

1. 326 

0.060 

0.220 

2.336 

2.  297 

2.203 

0.060 

0.280 

2.255 

2.214 

2.135 

0.070 

0.340 

2.160 

2.113 

2.038 

0.080 

0.440 

2.014 

1.963 

1.869 

0.110 

0.650 

1.709 

1.650 

1.535 

0.110 

0.700 

1.649 

1.591 

1.472 

0.170 

0.850 

1.383 

1.326 

1.236 

0.190 

0.975 

1.  184 

1.166 

1.165 

0.250 

0. 140 

2.124 

2.054 

2.101 

0.230 

0. 190 

2.100 

■ 2.021 

2.041 

0.250 

0.240 

2.008 

1.921 

1.953 

0.250 

0.310 

1.924 

1.830 

1 . 809 

0.260 

0.630 

1.523 

1.430 

1.388 

0.330 

0.800 

1.225 

1 . 145 

1.170 

0. 340 

0.090 

2.029 

1.971 

2.077 

I 
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Table  I continued. 


Composition  Bandgap  (eV) 


X 

y 

Eq  (3) 

Eq  (5) 

Measured  115] 

0.370 

0.120 

1.954 

1.883 

1.978 

0.410 

0.170 

1.843 

1.757 

1.850 

0.400 

0.030 

1.9  83 

1.957 

2.086 

0.460 

0.040 

1.887 

1.854 

2.019 

0.500 

0.040 

1.830 

1.797 

1.953 

0.570 

0.060 

1.723 

1.678 

1.854 

0.650 

0.060 

1.625 

1.580 

1.750 

0.710 

0.060 

1.558 

1.513 

1.686 

0.830 

0.100 

1.401 

1.346 

1.470 

0.870 

0.130 

1.333 

1.276 

1.389 

0.920 

0.190 

1.216 

1.166 

1 . 268 

Average 

error  (eV) 

-0.052 

-0,008 

rms  error  (eV) 

0.109 

0.118 

Calculation 

of  halt  ice  Con;;  t anj^ 

for  the  calculation  of  the  lattice  constant,  ttie 
throe  methods  discussed  above  are  identical,  as  given  in 
Hq.  (4).  for  the  calculation  of  the  bandgap , we  have 
elected  to  use  Moon's  pi'ocedure , since  it  has  some 
theoretical  basis,  whereas  the  other  two  are  more  or  less 
ad  hoc. 

The  qviaternary  lattice  constant  and  energy  bandgap 
C'.'ntours  ar'e  presented  in  figure  l(.i)-(i).  These  con- 
tours were  obtained  by  numi>rical  ^iolu^ions  ot  I’qs . (4) 
and  (b),  using  the  data  givi>n  in  Table  2.  In  all  ca'-.es, 
the  lowest  quaternary  bandgap  is  plotted  in  figure  1. 

The  shaded  regions  represent  compositions  for  which  the 
quaternary  alloy  is  an  itidirect  bandgap  material. 
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Estimation  of  Unknown  Ternary  Bowing  Constants 

Relatively  few  of  the  IIl-V  ternary  bowing  constants 
are  known  with  any  certainty,  especially  for  the  X(IOO) 
and  L(lll)  valleys,  figure  2 shows  the  reported  bowing 
parameters  for  several  ternary  materials  as  a function 
of  the  lattice  constant  difference  between  the  two  end- 
point binary  compounds.  While  there  is  some  uncertainty 
in  the  experimental  data  there  does  appear  to  be  a 
definite  trend  toward  larger  bowing  parameters  with 
larger  lattice  constant  differences.  The  solid  line 
relationship  shown  in  Figure  2 has  been  used  to  estimate 
several  unknown  bowing  constants,  as  identified  by  the 
asterisks  in  Table  2.  Since  bowing  is  known  to  occur, 
it  was  felt  that  this  procedure  is  better  than  arbitrarily 
setting  the  unknown  bowing  constants  to  zero. 


Summary 

Calculated  bandgap  and  lattice  constant  contours  have 
been  presented  for  nine  quaternary  III-V  material  systems. 
The  calculated  values  are  based  upon  interpolation  tech- 
niques which  use  known  values  of  ternary  III-V  parameters 
to  estimate  the  bandgap  and  lattice  constant  parameters  for 
the  quaternary  systems.  Since  there  is  some  uncertainty  in 
both  the  ternary  parameters  and  the  best  interpolation 
technique  these  calculated  values  must  be  considered  as 
first  order  approximations  until  more  experimental  data 
is  obtained  for  the  quaternary  systems.  However,  the 
curves  should  prove  useful  for  many  studies  of  the  general 
properties  of  the  increasingly  important  III-V  quaternary 
materials . 
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figure  1.  Energy  baiidgap  and  lattice  constant  contours  for  III-V  quaternary  alloys.  The 
solid  curves  are  the  energy  gap  contours  and  the  dashed  curves  are  lattice  con- 
stant contours.  The  shaded  region  shows  the  compositional  range  over  which  the 
material  has  an  indirect  bandgap. 
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Energy  bandgap  and  lattice  constant  contours  are  presented 
for  the  six  quaternary  alloys  formed  from  Al,  Ga,  In  and 
P,  As,  Sb,  with  compositions  of  the  form  A B Cj^_  _ D or 
AB  C Df  X • The  quaternary  handgaps  and  Ja^tice  constant: 
were^obtaiJied  using  an  interpolation  formula  proposed  by 
the  present  authors. 


Key  words;  quaternary  alloys,  bandgap,  lattice  constant 


Introduction 


In  a previous  paper  [1],  bandgap  and  lattice  constant 
contours  for  nine  III-V  quaternary  alloys  have  been  pre- 
sented. The  present  paper  extends  that  work  to  include  si: 
additional  quaternary  alloys.  Host  of  the  work  done  on 
III-V  quaternary  materials  has  been  on  one  or  more  of  the 
nine  quaternary  alloy  systems  consisting  of  two  elements 
from  group  III  and  two  elements  from  group  V.  These  are 
hereafter  referred  to  as  2:2  quaternary  alloys.  There  has 
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been  some  interest  [2,3,4]  in  the  set  of  IIl-V  qu<rternary 
alloy  systems  that  are  foi’roed  by  selecting,  either  three 
group  III  elements  and  one  group  V element,  or  one  group 
111  element  and  three  group  V elements.  These  are  here- 
after refeiTed  to  as  3:1  or  1:3  quaternary  alloys.  This 
paper  presents  energy  bandgap  and  lattice  constant  contours 
for  the  six  3:1  and  1:3  quaternary  alloys  formed  fix'm  Al, 

Ga , In  and  P,  As,  Sb.  The  calculations  are  based  on  an 
interpolation  procedure  which  is  similar  to  the  one  previ- 
ously applied  to  the  nine  2:2  quaternary  alloys  [1]. 


Not  at  Ion 


There  has  been  very  little  published  material  on  the 
3:1  and  1:3  III-V  quaternary  alloy  systems,  and  there  is 
no  generally-accepted  scheme  for  symbolically  describing 
the  composition  of  .such  material.s.  In  tliis  work  the 
following  convention  is  adopted:  the  composition  v.tri.ibles 
associated  with  the  three  variable  elements  are  in  the 
order  x,  y,  and  z.  These  are  constrained  by 

x+ytz^l.  (1] 

A 1:3  quaternary  alloy  will  be  described  as  AB  C D , whet'e 
A is  a group  111  element  and  B,  C and  P .ire  givu^  ^ 
elements  listed  accoi'iiing  to  increasing  atomic  number.  A 
3:1  quaternary  alloy  will  be  described  as  A^B  where 

A,  B and  C are  group  III  elements  listed  in  increasing 
atomic  number  and  D is  a group  V element.  A 3:1  or  1:3 
quaternary  alloy  par.imeter  is  described  by  a surface 
Q(x,y,z)  over  the  triangtilar  composition  pi  me  tO-ot*-!, 
0<y<l,  0<z<l;  xty+z=l)  shown  in  figure  1.  At  the  corners, 
the  values  of  the  parameter  for  the  three  biiuu’y  compounds 
Q(1,0,0)::B^,  Q(0,1,0)  ::  B^  and  Q(0,0,1)  = B^  are  obtained. 
Along  the  boumlaries  of  the  plane  the  paranieters  for  the 
three  ternary  alloys  Q(x,y,0)  = T^.,(u),  Q(0,y,z)  = 
and  Q(x,0,z)  ■ obtained,  where 

u = (l-xi-y)/2, 

V :=  (1-ytz)/?,  (2) 


w = (1  xtz)/2. 
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Fig.  1.  Illustrating  the  definitions  of  the 
ternary  (u,v,w)  and  quaternary  (x,y,2)  alloy 
composition  variables. 

Fstimation  oj^  A 1 Ij^  Farameters 

As  in  the  previous  paper  on  2:2  quaternary  alloys  ll], 
the  lattice  constants  of  the  teriuary  alloys  are  linearly 
interpolated  from  the  constituent  binary  conpcnents  as, 

T.  .(x)  = xB.  + (l-x)B.  , (3) 

1]  3 1 

and  the  bandgaps  of  the  ternary  alloys  are  obtained  using  a 
nonlinear  interpolation 

T..(x)  = xB.  + {l-x)B.  t C..x(l-x),  (*<) 

t]  3 t ID 

where  C..  is  the  bowing  parameter  for  the  ternary  alloy 
bandgap'^^' . . . The  quaternary  alloy  parameter  Q(x,y,2)  is 
interpolated  from  the  ternary  parameters  as 


mmmmmmmmmmrn 
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Q(x,y,z) 


^ * xzTj^gCw) 

xy+y z+xz 


(S) 


Here  the  T,,  are  the  ternary  alloy  bandgaps  or  lattice 
constants  given  by  Eq.  (3)  or  , respectively,  and  the 
variables  u,  v,  and  w are  as  defined  in  Eq.  (2),  The 
quaternary  parameters  calculated  according  to  Eq.  (5) 
reduce  to  the  appropriate  ternary  parameters  on  the  tri- 
angle boundaries  (x  or  y or  z equal  to  zero)  and  to  the 
average  of  the  three  ternary  parameters  at  the  center 
{x-y=z=l/3)  of  the  composition  plane. 


Calculation  of 

The  3:1  and  1:3  quaternary  lattice  constant  and  energy 
bandgap  contours  are  presented  in  Figures  2(a)-(f).  The 
contours  were  obtained  by  numerical  solution  of  Eqs . (3), 
(4)  and  (5),  using  the  data  given  in  Table  1 of  [1].  In 
all  cases,  tne  lowest  quaternary  bandgap  is  plotted  in 
Figure  2.  The  shaded  regions  represent  compositions  for 
which  the  quaternary  alloy  is  an  indirect  bandgap  material. 


Summary 

Energy  bandgap  and  lattice  constant  contours  for  the 
six  III-V  3:1  and  1:3  quaternary  alloy  systems  have  been 
presented.  Since  little  experimental  work  has  been  done 
with  these  materials,  verification  of  the  calculations 
cannot  be  made  at  the  present  time.  However,  the  calcu- 
lated parameters  should  be  reasonably  good  first  approxi- 
mations which  should  prove  useful  in  further  investiga- 
tions, realizing  that  Vegard's  law  seems  to  apply  for 
these  materials  C'*],  that  the  interpolation  procedure 
seems  to  be  accurate  for  the  2:2  quaternaries  [1],  and 
that  the  values  along  the  triangular  composition  boundary 
are  based  on  ternary  experimental  d.tta. 
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